April 2008 



Mirror dark matter and the new DAM A/LIBRA results 
A simple explanation for a beautiful experiment 



R. Foot 1 

School of Physics, 
University of Melbourne, 
Victoria 3010 Australia 



Recently, the DAMA/LIBRA experiment has convincingly confirmed the DAMA/Nal 
annual modulation signal, experimentally demonstrating the existence of non-baryonic 
dark matter in the halo of our galaxy. Meanwhile, in another part of town, other 
experiments such as CDMS and XENON10 have not detected any evidence for dark 
matter. One promising dark matter candidate which can reconcile the positive DAMA 
annual modulation signal with the null results from the other experiments, is mirror dark 
matter. We re-analyse the mirror matter interpretation of the DAMA annual modulation 
signal utilizing a) the new data from DAMA/LIBRA, including the measured energy 
dependence of the annual modulation signal b) an updated quenching factor which 
takes into account the channeling effect in Nal crystals and c) the latest constraints 
from CDMS/Ge, CDMS/Si and XENON10 experiments. We show that the simplest 
possibility of a He' (and/or H') dominated halo with a small O' component is sufficient 
to fully explain all of the dark matter experiments. We also point out that a certain 
class of hidden sector dark matter models, although theoretically less appealing and less 
constrained, can mimic the success of the mirror dark matter model and hence are also 
viable. 
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1 Introduction 



Recently the DAMA/LIBRA experiment has confirmedfl] the impressive DAMA/Nal 
annual modulation signal at around 8a C.L.[2]. Their signal is observed in the 2-6 keVee 
energy region, with periodicity and phase: 

T/year = 0.998 ± 0.003, t /day = 144 ±8, (1) 

which is beautifully consistent with the dark matter expectation of T/year = 1.0, 
to/day = 152.5. There are no known systematic effects which could produce the mod- 
ulation of the signal seen. Indeed, if there were some hypothetical systematic effect 
producing the modulation, it would be surprising if it had the same periodicity AND 
phase as the dark matter signal. Therefore, it is reasonable to believe that these exper- 
iments have detected galactic dark matter particles. 

Meanwhile other experiments have failed to detect positive evidence for dark mat- 
ter. Among the most sensitive of these, are the recent CDMS/Ge[3], CDMS/Si[4] and 
XENON10[5] experiments. Combined analysis of these null results with the positive 
results of DAMA appear to exclude standard WIMP explanations (such as the standard 
neutralino models) for the DAMA signal. Several alternative explanations of the DAMA 
signal have been proposed, such as the elastic scattering of light WIMPs (with mass ~ 
7 GeV)[6], and very light candidates (e.g. with mass in the keV range) scattering off 
electrons [7, 8]. While such phenomological models might be possible, it is a theoretical 
challenge to construct simple renormalizable extensions to the standard model which 
can accommodate these light particles. 

On the other hand, previous work[9, 10, 11, 12] has shown that mirror dark matter 
can explain the DAMA/Nal annual modulation signal, consistently with the null results 
of the other experiments. Mirror dark matter arises in simple renormalizable extensions 
of the standard model featuring a sector of particles and forces exactly isomorphic (or 
"mirror") to the known particles and forces. The simplest such model being the exact 
parity symmetric model[14] (for a review, see e.g. ref.[13]). These theories feature a 
spectrum of dark matter particles, A', with known masses. 

The success of the mirror dark matter theory in explaining the DAMA/Nal experi- 
ment consistently with the null results of the other dark matter detection experiments 
has to do with several of its key distinctive features: 

• Galactic halo mirror particles have a Maxwellian velocity distribution: f(u) = 
exp[ — u 2 / i>q] , however the velocity dispersion of heavy mirror particles {Ma> > 
M He ) is much less than the galactic rotational velocity: v (A') 2 <C v 2 ot . In standard 
WIMP models, v\ = v 2 rot . 

• The mirror dark matter mass spectrum consists of relatively light particles: Ma> < 
M Fe . Standard WIMP models typically feature M WIMP ~ 50 GeV. 

• Mirror dark matter interacts with ordinary matter via renormalizable photon- 
mirror photon kinetic mixing, leading to a recoil energy (E R ) dependent cross- 
section: oc -gr. Standard WIMP models have an approximately E R indepen- 
dent cross-section arising from a four-Fermion interaction. 
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These three features greatly increase the sensitivity of mirror dark matter experiments 
such as DAMA/Nal. The main reason is that the DAMA/Nal experiment is a lower 
threshold experiment than the sensitive CDMS and Xenon experiments, when the chan- 
neling effect is taken into account [15]. 

The purpose of this paper is to re-analyse the mirror matter interpretation of the 
DAMA annual modulation signal in the light of the new data from DAM A/LIBRA as 
well as the latest constraints from CDMS and XENON experiments. We also use an 
updated quenching factor, taking into account the channeling effect in Nal crystals, 
which was discussed recently by the DAMA collaboration [15]. 

The outline of this paper is as follows: in section II we give a brief review of mirror 
dark matter, which may be skipped by those already familiar with the subject. In 
section III, we give the necessary details relevant to direct detection experiments such 
as cross-section, halo distribution function etc. In section IV we give the detailed fit 
to the DAMA/LIBRA + DAMA/Nal combined results under the simple and plausible 
assumption that the heavy component of the galactic halo (i.e. with M > M He ) is 
dominated by just one mirror element. That is, one element in addition to the H' / He' 
predominant component. In section V we examine the constraints imposed by the 
negative results of the CDMS/Ge, CDMS/Si and XENON10 experiments. In section VI 
we point out that a class of hidden sector models can mimic the success of the mirror 
dark matter theory in explaining the experiments. Finally, in section VII we draw our 
conclusions. 



2 Mirror Dark Matter 

The exact parity symmetric model [14] is the minimal extension of the standard model 
which allows for an exact unbroken parity symmetry [x — > —x, t — > t]. According to this 
theory, each type of ordinary particle (electron, quark, photon etc) has a corresponding 
mirror partner (mirror electron, mirror quark, mirror photon etc), of the same mass. 
The two sets of particles form parallel sectors each with gauge symmetry G (where 
G = SU(3) <S> SU(2) <S> U(l) in the simplest case) so that the full gauge group is G ® 
G. The unbroken mirror symmetry maps x — > — x as well as ordinary particles into 
mirror particles. Exact unbroken time reversal symmetry also exists, with standard 
CPT identified as the product of exact T and exact P[14]. 

It has been argued that the stable mirror particles, mirror nucleons and mirror 
electrons are an interesting candidate for the inferred dark matter of the Universe (for a 
review, see Ref. [13]). Of course, to be a successful dark matter candidate, mirror matter 
needs to behave, macroscopically, differently to ordinary matter. In particular, four key 
distinctions need to be explained: 

• The cosmological abundance of mirror matter should be different to ordinary mat- 
ter, Qdark 7^ ^matter' 

• Mirror particles should give negligible contribution to the energy density at the 
epoch of big bang nucleosynthesis. 
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• Structure formation in the mirror sector must begin before ordinary matter radi- 
ation decoupling. 

• In spiral galaxies, the time scale for the collapse of mirror matter onto the disk 
must be much longer than that of ordinary matter. 

Clearly, mirror matter behaves differently to ordinary matter, at least macroscop- 
ically. It is hypothesised that this macroscopic asymmetry originates from effectively 
different initial conditions in the two sectors. The exact microscopic (Lagrangian) sym- 
metry between ordinary and mirror matter need never be broken. In particular, if or- 
dinary and mirror particles have different temperatures in the early Universe, T' <C T, 2 
then the mirror particles give negligible contribution to the energy density at the time 
of nucleosynthesis leading to standard big bang nucleosynthesis. Another consequence 
of T' <C T is that mirror photon decoupling occurs earlier than ordinary photon decou- 
pling - implying that mirror structure formation can begin before ordinary photon de- 
coupling. In this way, mirror matter-type dark matter can successfully explain the large 
scale structure formation (for detailed studies, see ref.[17, 18, 19]). Also, VL dark ^ f2 matter 
could also be due to different effective initial conditions in the early Universe (see e.g. 
ref.[20],[21] for some specific scenarios). 

If mirror matter is the inferred non-baryonic dark matter in the Universe, then the 
halo of our galaxy should be a gas of ionized mirror atoms and mirror electrons with 
a possible mirror star component 3 . Although dissipative, roughly spherical galactic 
mirror matter halo's can exist without collapsing provided that a heating mechanism 
exists - with ordinary supernova explosions being plausible candidates [24]. Obviously, 
the heating of the ordinary and mirror matter in spiral galaxies needs to be asymmetric, 
but again, due to different initial conditions in the early Universe, asymmetric heating 
is plausible. For example, the early Universe temperature asymmetry, T' <^ T (ex- 
pected from successful big bang nucleosynthesis and Large scale structure formation, as 
discussed above) implies that the primordial mirror helium/mirror hydrogen ratio will 
be much larger than the corresponding ordinary helium/ordinary hydrogen ratio [17]. 
Consequently the formation and evolution of stars in the ordinary and mirror sectors 
are completely different (see e.g. ref.[25] for some preliminary studies). The details of 
the evolution on (sub) galactic scales, is of course, very complex, and is yet to be fully 
understood. 

Ordinary and mirror particles interact with each other by gravity and via renormal- 
izable photon-mirror photon kinetic mixing: 4 

£ = \f" u f^ (2) 

2 It has been speculated that f«T might arise in certain inflation scenarios, see e.g. ref.[16]. 

3 The mirror star component can be probed by microlensing observations towards nearby galaxies [22]. 
The current situation is somewhat unclear, but a mass fraction of MACHO's around / <~ 0.1 is roughly 
compatible with the observations. For a review, see ref.[23] 

4 Technically, photon-mirror photon kinetic mixing arises from kinetic mixing of the U(l) and £7(1)' 
gauge fields, since only for the abelian U(l) gauge symmetry is such mixing gauge invariant [26]. The 
only other gauge invariant and renormalizable interactions mixing ordinary and mirror particles are the 
quartic Higgs - mirror Higgs interaction: Ac^V^'V and neutrino - mirror neutrino mass mixing[14, 27]. 
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where F^ v (F') is the field strength tensor for electromagnetism (mirror electromag- 
netism). One effect of photon-mirror photon kinetic mixing is to cause mirror charged 
particles (such as the mirror proton and mirror electron) to couple to ordinary photons 
with effective electric charge ee.[14, 28, 29] The various experimental implications of 
photon-mirror photon kinetic mixing have been reviewed in Ref. [30]. Of most relevance 
for this paper, is that this interaction enables mirror particles to elastically scatter off 
ordinary particles - essentially Rutherford scattering. 

A detector on Earth can therefore be used to detect halo mirror nuclei via elastic 
scattering. Several previous papers[9, 10, 11, 12] have explored this possibility, especially 
in view of the impressive dark matter signal from the DAMA/Nal experiment. The 
purpose of this paper is to re- analyse the mirror matter interpretation of the DAMA/Nal 
experiment using the new data from DAMA/LIBRA and also the latest constraints from 
the null results of the CDMS and XENON10 experiments. In the following section we 
give the technical details - the cross section and halo distribution function - necessary 
to facilitate this comparison of theory and experiment. 

3 Mirror dark matter implications for direct detec- 
tion experiments 

3.1 Cross-section and form factors 

Let us first briefly review the required technology (see ref. [9, 10, 11, 12] for more details). 
For definiteness, consider a halo mirror nuclei, A', of atomic number Z' scattering off 
an ordinary nucleus, A (in an ordinary matter detector) of atomic number Z. The 
cross-section is then just of the standard Rutherford form corresponding to a particle of 
electric charge Ze scattering off a particle of electric charge eZ'e. This cross-section can 
be expressed in terms of the recoil energy of the ordinary atom, E R , and the velocity of 
A' in the Earth's rest frame, v: 

^ X (3) 



dE R E\v 2 
where 

. 2ne 2 a 2 Z 2 Z' 2 
= W A F A^ r A)F 2 Aqr A i) (4) 

and Fx(qrx) {X = A, A') are the form factors which take into account the finite size 
of the nuclei and mirror nuclei, [q = (2MaEr) 1 I 2 is the momentum transfer and rx is 
the effective nuclear radius] 5 . A simple analytic expression for the form factor, which 
we adopt in our numerical work, is the one given by Helm[31, 32]: 

Fx(qrx) = 3 J -I^e~^ (5) 
qr x 



3 We use natural units, h = c = 1 throughout. 



4 



with r x = 1.14X 1 / 3 fm, s = 0.9 fm. In this equation, ji is the spherical Bessel function 
of index 1. 



3.2 Interaction rate and halo distribution function 

In an experiment such as DAMA/Nal, the measured quantity is the recoil energy, Er, 
of a target atom. The interaction rate is 

__ = ^N T n A , / — \v\drv 

dE R ^ J dE R k 

= }^ N T n A'-^ / rj— ; cru (6) 

where N T is the number of target atoms per kg of detector 6 . Also, n A > is the halo number 
density (at the Earth's location) of the mirror element, A' and /a'( v > v -e) A i s its velocity 
distribution {k = (tt^q) 3 / 2 is the normalization factor). Here, v is the velocity of the 
halo particles relative to the Earth, and is the Earth's velocity relative to the galactic 
halo. [Note that the bold font is used to indicate that the quantities are vectors]. The 
lower velocity limit, v m i n (En), is given by the kinematic relation: 



(Ma + Ma') 2 E r 
\ 2M A M% ' 1 j 

Considering a particular mirror chemical element, A' (e.g. A 1 = H', He', O' etc), the 
velocity distribution for these halo mirror particles is then: 



fA'(v,v E ) = exp 



•^M A ,(v + v s ) 2 /T 



= exp[-(v + v E ) 2 /v%] (8) 



where v 2 = 2T/M A ,. 

The assumption of approximate hydrostatic equilibrium for the halo particles implies 
a relation between T and the local rotational velocity, f rot :[10] 



fj,M p v, 



t = ^y ot (9) 

where fiM p is the mean mass of the particles comprising the mirror (gas) component of 
the halo (M p is the proton mass). Note that the Maxwellian distribution should be an 
excellent approximation in the case of mirror dark matter, since the self interactions of 
the particles ensure that the halo is locally thermalized. 
The velocity integral in Eq.(6), 

i(Er) s r ^^<?v no) 

J\v\>v min {E R ) k\V\ 



6 For detectors with more than one target element we must work out the event rate for each element 
separately and add them up to get the total event rate. 
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is standard (similar integrals occur in the usual WIMP interpretation 7 ) and can easily be 
evaluated in terms of error functions assuming a Maxwellian dark matter distribution [32], 

fA'(v,V E )/k = M)" 3 / 2 exp[-(y + V E ) 2 /^\, 

I(Er) = ^ [erf(x + y) - erf(x - y)\ (11) 

where 

xe^!,,^. (12) 

^0 ^0 

The Earth's velocity relative to the galactic halo, ve, has an estimated mean value of 
(v E ) — v rot + 12 km/s, with v rot , the local rotational velocity, in the 90% C.L. range[33], 

170 km/s ~ v rot ~ 270 km/s. (13) 

While some estimates put narrower limits on the local rotational velocity, it is useful 
to allow for a broad range for v rot since it can also approximate the effect of bulk halo 
rotation. 

As can be seen from Eq.(8,9), in the case of mirror matter-type dark matter, the Vq 
value for a particular halo component element, A', depends on the chemical composition 
of the halo. In general, 

V 2 rot M A , ' 

The most abundant mirror elements are expected to be H', He', generated in the early 
Universe from mirror big bang nucleosynthesis (BBN). Heavier mirror elements are 
expected to be generated in mirror stars, or possibly, in the early universe if mirror BBN 
occurs early enough, so that the number densities are large enough, to allow three-body 
processes such as the triple alpha' process to be effective. It is useful, therefore, to 
consider two limiting cases: first that the halo is dominated by He' and the second is 
that the halo is dominated by H' . The mean mass of the particles in the halo are then 
(taking into account that the light halo mirror atoms should be fully ionized): 

fiM p ~ Mflv/2 ~ 0.5 GeV for H' dominated halo, 

fiM p ~ M He ,/3 ~ 1.3 GeV for He' dominated halo. (15) 

The v values can then easily be obtained from Eq.(14): 



Vo(A') = v (H')J ^j- ~ ~7=\ TF~~ km/s for H' dominated halo 
V Ma> v2 V Ma' 

v (A') = v (He')\ r^f e -ry^- km/s for He' dominated halo. (16) 



7 However in the WIMP case the upper velocity limit is finite, corresponding to the galactic escape 
velocity. While for mirror dark matter, the upper limit is infinite due to the self interactions of the 
mirror particles. 
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Mirror BBN[17] suggests that He' dominates over H' . It happens that the recoil energy 
thresholds of the DAMA/Nal and the other experiments are sufficiently high that these 
experiments are only directly sensitive to mirror elements heavier than about carbon, 
which means that vq{A!) <C v rot for these elements - independently of whether He' or 
H' dominates the halo. It turns out that our main results (such as the allowed region 
in figure 4) do not depend very significantly on whether we assume that He' or H' 
dominates the mass of the Halo. 

Finally, we see from Eq.(6) that the interaction rate, dR/dEn, is proportional to 
e 2 riA', where e is the photon-mirror photon kinetic mixing parameter and ua' is the halo 
number density of the species A'. The dark matter halo mass density at the Earth's 
location is approximately 0.3 GeV/cm 3 , so that the halo mass fraction, £4', of the species 
A', is thus: £4/ = ha'Ma' /(0.3 GeV/cm 3 ). Hence, we find that the interaction rate 
is proportional to (ty/^) 2 - We call the quantity, ey^^, the cross-section abundance 
coefficient, which we can determine from the amplitude of the DAM A annual modulation 
signal. The only other variables in the model are the rotation velocity, v rot , which is 
expected to be in the range, Eq.(13), and the spectrum of heavy mirror elements, A'. 

4 Mirror matter interpretation of DAMA/Nal & 
DAMA/LIBRA annual modulation signal 

The DAMA experiments employ large mass scintillation detectors located in the Gran 
Sasso underground laboratory. The first generation, DAMA/Nal experiment ran for 
7 years, presenting final results in 2003 [2]. The second generation, DAMA/LIBRA 
experiment has been running since 2003, recently presenting its first results from 4 
years of data[l]. 

The DAMA/Nal and DAMA/LIBRA experiments are very sensitive to mirror matter- 
type dark matter because of the relatively low energy threshold of 2 keVee. The unit, 
keVee is the so-called electron equivalent energy, which is the energy of the event if due 
to an electron recoil. The actual nuclear recoil energy (in keV) is given by: keVee/q, 
where q is the quenching factor. 

The DAMA/Nal and DAMA/LIBRA experiments extract their signal by using the 
annual modulation signature [34], which arises because of the Earth's orbital motion. 
The point is that the interaction rate, Eq.(6), depends on ve, which varies due to the 
Earth's motion around the sun: 



where v® = v rot + 12 km/s ~ 230 km/s is the sun's velocity with respect to the galactic 
halo and v e ~ 30 km/s is the Earth's orbital velocity around the Sun (t = 152.5 days 
and uj = 2tt/T, with T = 1 year). The inclination of the Earth's orbital plane relative 
to the galactic plane is 7 ~ 60°, which implies that Ave — 15 km/s. The differential 
interaction rate, Eq.(6), can be expanded in a Taylor series around ve = u©, leading to 



v E (t) 



v Q + v® cos 7 cos u{t — t ) 
v Q + Ave cos uj(t — to) 



(17) 
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an annual modulation term: 



Ri = R° + R] cos u(t - t ) 



(18) 



where 



BS = 



R) 



I rEi+AE / 

-I 

\E Je 



dR 



dE R 

V E =VQ 

rEi+AE Q t dR\ 

ae .//;, ! ' 



dE 



R, 



dv E \dE R/ 



Av E dE R 



(19) 



ve=vq 



According to the DAMA analysis[l], they indeed find an annual modulation at more 
than 8a C.L. Their data fit gives T = (0.998 ± 0.003) years and t = 144 ± 8 days, 
consistent with the expected values. [The expected value for to is 152.5 days (2 June), 
where the Earth's velocity reaches a maximum with respect to the galactic halo] . Their 
signal occurs in the 2-6 keVee energy range, with amplitude 



R 1 = (0.0129 ± 0.0016) cpd/kg/keVee [cpd = counts per day]. 



(20) 



The theoretical rate for the modulated and absolute component can be evaluated in 
our model from Eqs.(19,6,10,ll), leading to: 



dE? 
dE R 

dR 1 



dE 



R 



E 

A' 

E 



N T n A '\I(E R ,y ) 



El 



N T riA'\A.y 
El 



'or 
9y, 



(21) 



y=yo 



where y = (v E )/v , Ay = Av E /v , and 

I(E R ,y ) 



'df\ 



yo 



+ 



V^voyo 



e -(x-y ) 2 _|_ e -(x+y) 2 



(22) 



Recall, from Eqs.(12,7,14) that x oc y/E^, while y ~ = J^jfc, so that y ~ 3.4 



M 



A' 



5.5 



Mr 



for a He' (H') dominated halo. 

Consider the interactions of one halo component A'. What is the behaviour of the 
annual modulation amplitude, dR l /dE R) as a function of E R 7 At low E R , (where 
i C I/), dR 1 /dE R is negative, as E R increases, dR 1 /dE R changes sign, and reaches a 
maximum when 



x(E R ) ~ yo. (23) 

At high E R (x > y), dR 1 /dE R -> 0. 

Of course, in order to do a detailed fit of the experimental data to a given theoretical 
model, we must take into account the finite resolution of the DAMA detector and the 
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physics of the quenching factor. In all of our numerical work both of these effects are 
included. The experimental energy resolution is taken into account by convolving the 
rate with a Gaussian, with a obtained from ref.[35]: cr/E = a/ \J E{keV) + f3 with 
a = 0.448 ± 0.035, (3 = (9.1 ± 5.1) x 10~ 3 . 

The quenching factor is more complicated. Our previous work assumed q^ a = 0.3 
and qi = 0.09. However, recently Drobyshevski[36] has pointed out the possible impor- 
tance of the 'channeling effect' in Nal crystals, which means that ions moving approx- 
imately in the direction of crystal axes and planes get 'channeled' whereby they give 
an unquenched signal, q ~ 1. A detailed study and modelling of this effect has been 
done recently by the DAMA collaboration [15]. According to their model, the fraction 
of events, fNa,i, with quenching factor q ~ 1, in the low recoil energy region of interest 
is approximately: 

fNa ~Y+\lAEdk^)' fl ~ 1 + 0.75E R (keV) ^ 

The remaining fraction of events have quenching factors, q^ a ~ 0.3, qj ~ 0.09 (but 
are broadened by straggling). The channeling effect is very important, as it effectively 
reduces the threshold of the DAMA experiments down to an actual recoil energy of 2 
keV. It turns out that the fitted cross-section is over an order of magnitude smaller 
compared to the case where the channeling effect is ignored. This means that the 
interaction rate is dominated by the channeled events. 

The 2 keV threshold is still too high for the DAMA experiment to be very sensitive 
to the dominant, He' or H' component, except possibly for effects in the E R ~ 2.5 keV 
region from the tail of the He' induced events. Basically, these nuclei are too light to 
give much of a signal above the DAMA/Nal energy threshold. [Of course, the DAMA 
experiment is indirectly sensitive to the H' /He' component, which as we discussed in 
section 3, has the important effect of reducing the value of vq(A') ]. DAMA is only 
directly sensitive to mirror nuclei heavier than about carbon. In previous work we 
considered having several components, such as O' , Si' , Fe' in ref.[12], however it turns 
out that this is not necessary to fit the new DAM A/LIBRA data. The reason is two 
fold: first, DAMA/LIBRA have found a lower annual modulation amplitude in the 4-6 
keVee energy region than was found in the earlier DAMA/Nal experiment, second, the 
channeling effect also modifies the predicted spectrum. The net effect, as we will see, is 
that the interaction of a single halo element predicts a signal of the correct shape in the 
2-6 keVee range. We therefore assume this simplest possibility, that is, that the spectrum 
of heavy elements (in this context 'heavy' means heavier than M# e ) is dominated by 
just one element, A'. This would also approximate a narrow range of elements. The 
detailed predictions will depend on the mass of the element A', the rotational velocity 
v rot and the cross-section abundance coefficient, ey^T. 

Note that the cross-section abundance coefficient, {e^/CA') 2 , is just an overall normal- 
ization factor, and doesn't change the shape of the annual modulation energy spectrum. 
This factor can be adjusted so that R 1 fits the measured amplitude averaged over the 
2-6 keVee energy range. Of course, only the amplitude is a free parameter, the phase, 
periodicity and the approximate cosine form of the signal are all predicted, as they are 
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for many other dark matter models. The fit of the annual modulation data is shown in 
figure 1. 
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Figure 1: Measured single hit residuals for the (2-6) keVee region (obtained from figure 
6 or ref.[l]) compared to the predicted cosine modulation with phase t /day = 152.5 
day and T/year = 1. 

An important new result from DAM A/LIBRA is that we now have a much better 
idea of the shape of the annual modulation energy spectrum. The DAMA collaboration 
separate the data in the 2-6 keVee range into 8 energy bins and fit the modulation in 
each energy bin to the predicted cosine distribution, obtaining modulation amplitudes, 
S™ ± <7j, for each energy bin. These amplitudes can be compared to the theoretically 
computed quantities, R}, taking into account the experimental resolution and quenching 
factors as described earlier. We can compute a x 2 quantity, 

X 2 = E(^ - Rl) 2 /*- (25) 
i=l 

which is a function of the three parameters: v ro t, ^\f^A' and the identity of the dominant 
element A'. Thus, with eight independent data bins we have five degrees of freedom. We 
find numerically that for each value of v rot [in the range, Eq.(13)], the model provides a 
good fit to the data, with: 

Xmin ~ 5.5, for H' dominated halo 

Xmin ~ 3.5, for He' dominated halo (26) 
10 




for five degrees of freedom. The data does constrain the other two parameters. Numer- 
ically we obtain the 3<r range for ey/^: 

\ e \J^L ~ (1.0±0.3)xl0- 9 f V \ ot - V (27) 
11 V 0.1 V ; \210 km/s J V ; 

for both H' and He' dominated halo. 

The allowed range for A' also depends on v rot . We find that we can fit the DAMA 
data for A' ranging in mass from oxygen to sulfur, i.e. 15 GeV ~ M A > ~ 30 GeV. In 
figure 2 we have given the best fits for four illustrative examples, A' = 0', v rot = 280 
km/s, A' = Ne\ v rot = 230 km/s, A' = Mg', v rot = 195 km/s and A' = Si', v rot = 170 
km/s. Figure 2a assumes a H' mass dominated halo, while figure 2b assumes a He' 
mass dominated halo. 
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Figure 2a: Energy distribution of the cosine modulation amplitude for four illustrative 
cases: A' = Si',v rat = 170 km/s (solid line) A' = Mg',v rot = 195 km/s (long-dashed 
line), A' = Ne',v rot = 230 km/s (short-dashed line), A' = 0',v rot = 280 km/s (dotted 
line). In each case e^/^ is fixed so that the mean amplitude is 0.0129 cpd/kg/keVee. 
Also shown are the DAMA/Nal & DAMA/LIBRA combined data from figure 9 of ref. [1]. 
This figure assumes that the mass of the halo is dominated by H'. 
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Figure 2b: Same as figure 2a, except the mass of the halo is assumed to be dominated 
by He'. 

The shape of the measured annual modulation energy spectrum agrees very well 
with the mirror dark matter theory. It is essentially a two parameter fit, since as figure 
2 demonstrates, there is an approximate redundancy. Once the height and position of 
the peak is fitted, there is little freedom left to change the shape of the distribution. 
Thus, the fact that the measured shape agrees with that predicted by the mirror matter 
theory is a significant test of the theory. We illustrate this in figure 3 by fixing v rot and 
varying around the best fit value, for a particular example. 
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Figure 3: Annual modulation amplitude for the case of A' = O', (solid line), A' = C 
(dashed line), A' = Ne' (short-dashed line), ie. with A' — 16 ± 4. In each case we have 
used the same v rot = 280 km/s, assumed a He' mass dominated halo, and normalized the 
height of the peak to be the same. This figure illustrates that the position of the peak 
determines A', with the shape of the energy spectrum being approximately unchanged. 



As figure 3 illustrates, the position of the peak provides a measurement of the mass of 
the element A'. Allowing for the possible uncertainty in v rot , pins down the mass range 
for the element A' to the range: 15 GeV ~ M^> ~ 30 GeV. It is interesting that this 
mass range is theoretically consistent [Mue < Ma> < Mp e ) and also close to the naive 
expectation, of oxygen mass. Oxygen is the most abundant heavy ordinary element in 
the Universe, and it would not be surprising if the mirror sector were similar. 

If we consider A' as a parameter, then we can evaluate the 3a (5a) DAMA allowed 
region by finding the contours in (v rot , A') space where x 2 = Xmin + 9 (x 2 = Xmin + 2 5), 
where at each point, we vary ey^T such that the normalization is correct. We plot this 
DAMA allowed region in figure 4, for the case of a He' mass dominated halo. The case 
of a H' dominated halo is very similiar. 
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Figure 4: The 3<r (dashed line) and 5<r (solid line) allowed region of v ro t, A' parameter 
space consistent with the measured energy distribution of the modulation amplitude. 
In each case we have fixed e^/^ so that the mean modulation amplitude is 0.0129 
cpd/kg/keVee. 



The main features of our numerical results can be approximately understood: The 
data shows evidence for a peak in the annual modulation energy spectrum. A peak 
is expected analytically. As we discussed earlier [Eq.(23)], the position of the peak 
corresponds to the value of Er (which we denote as E p R eak ) where x ~ yo, or equivalently, 
to the value of Er where v r , 
we find: 



(E R eak ) ~ v rot . Evaluating this expression, using Eq.(7), 



M A , 



M A 



Vrot 



2M A 



(28) 



This shows again that the value of E R at the maximum of dR 1 /dE R , provides a deter- 
mination of the mass of the dominant heavy mirror element A'. The experimental data 
suggests that the maximum occurs roughly at Er ~ 2.75 keVee. Evaluating Eq.(28) for 
A — I and for A = Na gives: 



M A , 
M p 

M A , 
M p 



13.7 



Vrot 



( ;«)(; 

5.8 



) 



(300fem/s) °' 25 



, for A = I 

Na 



0.107 

, for A 



(29) 



Thus, we see that our allowed region corresponds to the case where A' is interacting 
predominately with I. We might have expected a second distinct region of allowed 
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parameter space - a second branch at lower v rot - where the rate is dominated by in- 
teractions with Na. Numerically it turns out that the Na peak is washed out by the 
interactions of /, whose peak then occurs around 1.1 keVee (if the peak in Na is at 2.75 
keVee), and is smeared out into the E > 2 keVee range by the detector resolution. The 
height of the / peak is also enhanced c.f. with Na peak, because a) the basic cross- 
section on / is larger than Na because of the larger / nuclear electric charge, b) the 
I peak occurs at lower Er and oc -^r, c) the fraction of channeled events strongly 

increases towards unity as Er — > in the adopted quenching model of ref.[15]. 

Finally, in figure 5, we give the absolute (unmodulated) spectrum, predicted by 
the model, and compare it to the DAMA spectrum. Of course, the measured absolute 
spectrum has a significant background component which is removed when the time 
varying component is extracted. Figure 5 indicates that in this model the required 
background spectrum is roughly consistent - it is smooth and positive. There is also a 
hint of the large rise in absolute rate predicted by the model at low energies. Of course, 
this is below the experimental threshold, and this hint should not be taken too seriously. 

Note that because of the sharply rising absolute spectrum at low E R our predicted 
shape of the absolute spectrum is somewhat sensitive to uncertainties in detector reso- 
lution and quenching factor. The shape of the modulation spectrum, on the other hand, 
is much less sensitive to these uncertainties. 
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Figure 5a: Predicted absolute (unmodulated) rate, j^, for the DAMA/Nal experiment. 
We take the same four examples of figure 2a: A' = Si',v ro t = 170 km/s (solid line) 
A' = Mg',v rot = 195 km/s (long-dashed line), A' = Ne',v rot = 230 km/s (short-dashed 
line), A' = 0',v rot = 280 km/s (dotted line). As with figure 2a, these results assume a 
H' dominated halo by mass. Also shown is the experimental (single hit) data. 
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Figure 5b: Same as figure 5a, except the mass of the halo is assumed to be dominated 
by He'. 



The value of e^/£j> obtained from the fit to the DAMA/Nal and DAM A/LIBRA 
experiments suggests a value for e in the range 10~ 9 — 10~ 8 for 10~ 3 < £a' < 10 _1 . 
Values of e in this range have many interesting applications, see e.g. ref. [37, 38]. A 
mirror sector, interacting with the ordinary sector, with kinetic mixing in this range is 
also consistent with Laboratory [39] and big bang nucleosynthesis constraints [40]. 

5 Constraints from the null experiments: CDMS / Ge, 
CDMS/Si, XENON10 

The annual modulation signal has only been seen in the DAMA/Nal and DAMA/LIBRA 
experiments. Other direct detection experiments have not found any dark matter signal. 
The most sensitive of these experiments are the CDMS and XENON10 experiments. 
These experiments use sensitive background elimination techniques and place limits on 
the absolute dark matter event rate, above their recoil energy detection thresholds. 

Importantly, these experiments are all higher threshold experiments, with claimed 
nuclear recoil energy thresholds of 10 keV for CDMS/Ge, 7 keV for CDMS/Si and 4.5 
keV for XENON10. The fact that mirror dark matter is light {Ma* ~ Mo) has a narrow 
velocity dispersion, vo{A') 2 <C vf ot , and has cross-section proportional to l/E^ all help 
to enhance the signal in DAMA/Nal relative to these higher threshold experiments. 

We have evaluated the contraints on our model from the CDMS/Ge[3], CDMS/Si[4] 
and XENON10[5] experiments. The published net effective exposures for these experi- 
ments are 121 kg-d, 12 kg-d and ~ 100 kg-d for CDMS/Ge, CDMS/Si and XENON 
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respectively. As, indicated by figure 5, we expect a sharply falling absolute event rate. 
Therefore the most sensitive energy region in these higher threshold experiments will 
be near their energy threshold. In the threshold region, CDMS/Ge[3], CDMS/Si[4] and 
XENON10[5] found 0, and 1 event respectively. The Poisson 95% exclusion region is 
then N = 3 for the CDMS experiments and N = 4.74 for the XENON experiment. We 
have computed these exclusion limits for the parameter space of interest for the DAMA 
experiment, which we show in figure 6. As the figure clearly demonstrates, our favoured 
region of parameter space is consistent with the null results of these experiments. 
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Figure 6: DAMA allowed region, as shown in figure 4, except with the addition of the 
95% C.L. exclusion limits from the XENON10, CDMS/Ge and CDMS/Si experiments 
(the regions above the exclusion contours are the disfavoured region). 

Finally, note that there are a few experiments with lower thresholds than DAMA/Nal. 
In particular there is the CRESSTI experiment with threshold of 0.6 keV[41]. The com- 
patibility of the CRESST experiment with the DAMA/Nal annual modulation signal 
was analysed in the mirror matter context in ref.[9], but without taking into account 
the channeling effect in the DAMA analysis. It was found that the predicted rate for 
the CRESST experiment was about a factor of two higher than the measured rate. In- 
cluding the channeling effect in the DAMA analysis lowers the predicted cross section 
by more than an order of magnitude, which means that the predicted event rate for 
the CRESST experiment is around 10% of their measured event rate. Thus the mir- 
ror dark matter interpretation of the DAMA/Nal signal is completely consistent with 
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the results of the low threshold CRESST experiment. Another low threshold experi- 
ment of note is the TEXONO experiment [42] with a claimed threshold of 0.2 keV. The 
low exposure and large background currently limit the usefulness of that experiment. 
However, the sharply rising event rate at low recoil energies, predicted by the mirror 
dark matter explanation of the DAMA/Nal experiment, might potentially be probed 
by future upgrades of experiments such as TEXONO and CRESST. The possibility of 
searching for annual modulation in such low threshold experiments would, of course, 
also be worthwhile. 

Finally note that there are earlier CDMS/Ge experiments with thresholds of 5 
keV[43] and 7 keV[4]. However the lower exposure and high backgrounds of these 
early runs make these experiments less sensitive than the XENON10 experiment for 
the mirror dark matter interpretation of the DAM A experiment, and therefore, within 
this theoretical framework, they do not disfavour any of the DAMA allowed parameter 
space. 

6 Hidden sector dark matter models 

Mirror matter arises if one assumes an extra sector of particles and forces exactly isomor- 
phic to the known particles. With the assumption that T' <C T in the early Universe, 
this simple model is compatible with the inferred properties of dark matter, such as its 
large scale structure, and, as we have shown, with all of the direct detection experiments 
as well. 

Of course, one can imagine generic hidden sector models which can approximately 
mimic the effects of the mirror matter model for dark matter experiments. The simplest 
such model would have a) two stable particles with masses M x and M 2 , such that 
M\ <C M 2 and b) interact with each other via an unbroken U(l)' (and possibly other 
gauge interactions), and c) interact with ordinary matter via U(l)y <S> U(l)' kinetic 
mixing. With these assumptions, one can have v (M 2 ) 2 <C v^ ot if the mass density of the 
halo is dominated by the lighter particles of mass Mi. Provided that 15 GeV ~ M 2 ~ 30 
GeV such a model can give similar predictions to the mirror matter theory, and thus is 
experimentally viable - at the moment. 

However it is interesting that all these features automatically occur in the mirror 
matter theory, which is theoretically tightly constrained. For example, in the mirror 
matter theory, M 2 is predicted to be in the range Mn e < Ma> < Mp e . Thus, the mirror 
dark matter theory is theoretically favoured over some generic alternative hidden sector 
model. This theoretical prejudice will be put to the test by future experiments, which 
we await with interest. 

7 Conclusion 

Mirror dark matter arises in simple and renormalizable extensions of the standard model 
whereby an exactly isomorphic sector of particles and forces is hypothesized which in- 
teract with the known particles via renormalizable U(l) kinetic mixing. Such extensions 
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of the standard model are theoretically well motivated, allowing space-time parity to 
exist as an exact unbroken symmetry. They also supply a spectrum of necessarily stable 
dark matter candidates of known masses. 

We have demonstrated that this dark matter theory can simply and fully explain the 
annual modulation signal seen in the DAMA/Nal and DAMA/LIBRA experiments. Our 
results indicate that the simplest and most plausible scenario: a predominant H' /He' 
halo with a small subdominant ~ O' component is sufficient to explain the data. The 
required photon-mirror photon mixing parameter is ~ 10~ 9 , which is consistent with all 
other experiments, and observations. In particular, we have shown that this simple sce- 
nario is also completely compatible with the latest null results of the other experiments, 
such as CDMS and XENON10. 

Generic hidden sector models which approximately mimic the effects of mirror matter 
type dark matter for direct detection experiments are also possible, but theoretically 
less appealing. Future experiments will, of course, be the final arbitrator. 
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